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Quantum phase transitions (QPTs) at zero tempera- 
ture are characterized by the change in the properties 
of the ground state of a many-body system caused by 
modifications in the interactions among its constituents 
Q. QPTs are induced as a parameter g in the system 
Hamiltonian H(g) is varied across a point g c . These 
phase transitions are completely driven by the quantum 
fluctuations and are incarnated via the non-analytic be- 
haviors of the ground-state properties at the transition 
points. On the other hand, as it is well known, the con- 
cept of entanglement lies at the heart of quantum me- 
chanics 0- Q Therefore, one expects that quantum en- 
tanglement (QE) should play an important role in QPTs. 
Recently, a great deal of effort has been devoted to un- 
derstanding their connection. Indeed, it has been ob- 
served that the quantum phase transitions are signaled 
by critical behaviors of the concurrence, a measurement 
of bipartite entanglement [4|, in a number of spin models 
llHSEIllIilillllEQIli However, they are 
not universal. For example, for the transverse-field Ising 
model, Osterloh et. al. found that the first order deriva- 
tive of the concurrence diverges at the transition point 
and obeys a scaling law in its vicinity (jj 0. On the 
other hand, for the antiferromagnetic XX Z chain (ID), 
the concurrence behaves in a completely different way. It 
is a continuous function of the anisotropic parameter and 
reaches its maximum at the transition point || . While in 
two- and three-dimensional (2D & 3D) XX Z models [T(ij. 
it develops a cusp- like behavior around the critical point. 
Another interesting example is the so-called J\~J 2 model 
(or Majumdar- Ghosh model) 0. Its ground state un- 
dergoes a first-order phase transition at J2/J1 = 0.5 and 
a continuous phase transition at J2/J1 — 0.241. How- 
ever, unlike the previous cases, the concurrence itself is 
now discontinuous at 0.5, while at 0.241, no discernable 
structure has been found |15| . Therefore, a natural ques- 
tion arose is why the same quantity, which measures en- 
tanglement between two localized spins, has such differ- 
ent behaviors, such as singularity [l3j. maximum, scal- 
ing, etc., at QPTs? Further, is there a unified picture 
of QE at the QPT? Obviously, investigation on these is- 
sues will not only deepen our understanding in QPTs 
but also strengthen the connection between condensed 
matter physics and quantum information theory |17| . 

In this Letter, we study these issues by detail anal- 



ysis of several well known spin models aim at giving a 
unified picture of QE at QPTs. For definitcness, wc 
choose the concurrence, C, as the measure of pairwise 
entanglement Q in this work. We characterize the be- 
haviors of the concurrence into three types and empha- 
size the important role played by the low-lying excitation 
spectra reconstruction of many-body systems around the 
transition points in determining the critical behaviors 
of the concurrence. More precisely, we show that the 
low-lying excitation spectra of these models are recon- 
structed in three qualitatively different ways around g c 
and hence, their concurrences show the above-mentioned 
non-universal behaviors at the QPT points. Therefore, 
we are able to give a unified and intuitive picture to un- 
derstand the relation between the QPT and the entan- 
glement. 

To have concrete discussions, we concentrate on three 
well studied quantum spin models in one dimension, the 
,/i — J-2 model, the XX Z model, and the transverse-field 
Ising model, defined by the following Hamiltonians: 
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In these Hamiltonians, sf , and sf denote the spin oper- 
ators at lattice site i. J±, J 2l A and A > are interaction 
parameters. The periodic boundary condition si = sjv+i 
is assumed. In Table I, we present basic features of these 
three models. The concurrence of the J\ — J 2 model and 
the spin ladder is exemplified in Fig. ^ while the con- 
currence of the Ising model and the XX Z model can be 
found in Refs. 0,01 and Refs. pi I9I IT(i| respectively. As 
summarized in Table I, there exist three types of QE at 
QPT points: (I) C is discontinuous; (II) C is continuous 
and exhibits maximum at QPTs; (III) C is continuous at 
QPTs but its higher order derivative exhibits extremum. 
We shall show that these seemingly different critical be- 
haviors of the concurrence can be understood on the ba- 
sis of the low-lying spectra reconstruction of these models 
around their QPT points. For this purpose, we present 
the ground-state energy and some low excited-state en- 
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TABLE I: The basic features of typical spin models, such as the properties of the concurrence, level-crossing(LC) in the 
ground-state (GS) and the first excited-state(ES), symmetry at the transition point, and type of phase transition. 
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Model (QPT point) 


GS LC 


ES LC 


concurrence 


symmetry 


transition type 


type 


XXZ chain(A = -1) 


Yes 




singular 






I 


Ji — J2 model( J2 = 0.5) 


Yes 




singular 






I 


XXZ chain(A = 1) 


No 


Yes 


maximum, not singular 


SU(2) point 


order-to-order 


II 


spin ladder( J = 0) 


No 


Yes 


maximum, not singular 


SU(2)®SU(2) 


disorder-to-disorder 


II 


XXZ 2&3D(A = 1) 


No 


Yes 


maximum, singular 


SU(2) 


order-to-order 


II 


Ji - J 2 model( J 2 ~ 0.241) 


No 


Yes 


not maximum 


unknown 


order-to- disorder 


III 


Ising model(A = 1) 


No 


No 


singular, not maximum 


unknown 


disorder-to-order 


III 
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FIG. 1: Concurrence of a 8-site Ji — J2 model (LEFT) and 
a 2 x 4 spin-ladder system (RIGHT). 



FIG. 3: LEFT: The second derivative of the concurrence with 
respect to J2 for various-size systems. RIGHT: The scale of 
the minimum point of the derivatives of the concurrence with 
different order. 
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FIG. 2: Low-lying energy spectra of four typical spin models. 



ergies of the J\ — J2 model, the XXZ model, the spin 
ladder model and the transverse-field Ising model, on a 
finite chain in Fig. [21 

For type I, we see immediately that, in both the J\ — J2 
and the XXZ models, there exist ground-state level- 
crossing at J2/J1 = 0.5 and A = — 1, respectively. Since 
the concurrence is a measurement of bipartite entangle- 
ment in the ground state, it is not difficult to see why it 
changes discontinuously at these transition points. The 
same behavior was also observed in the first order tran- 
sition in the spin model with mutual exchanges [TT| . 



To differentiate type II and III, we must take low-lying 
excitations into account. Let us study the spectra of the 
XXZ model at A = 1, the spin ladder model at J± = 0, 
and the transverse-field Ising model at A = 1. It has been 
proven that the ground states of these models are non- 
degenerate around the corresponding transition points 
[lfjj . At the first glance, it seems that this observation 
does not help us very much in understanding the critical 
behaviors of the concurrence in these models. However, 
as it is well known, QPTs are not solely dictated by the 
ground state of a specific model, especially on a finite sys- 
tem. They also depend on interconnection between the 
ground state and the low-lying excited states of the sys- 
tem P. In other words, the change of the ground-state 
properties is greatly affected by matrix elements of the 
relevant collective-mode operators, such as the spin-wave 
operators or the particle-density- wave operators, which 
relates the low-lying excited states to the ground state. 
In particular, in Ref. j2jj, we pointed out that the QPT 
is actually induced by the reconstruction of the excita- 
tion energy spectra. As a consequence, when a transition 
takes place between two ordered phases of a system, such 
as the one for the XXZ model at A = 1, which separates 
the XY regime from the Ising regime, a level-crossing in 
its excited spectrum must occur. In fact, in Fig. [21 we 
see clearly the existence of such a level-crossing between 
the first and second excited states at A = 1 for the XXZ 
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model and at J± = for the spin ladder model. Both 
excited states have total spin 5=1 and are connected 
to the ground state, which is a spin singlet, via the spin- 
wave operators. It also implies that the system enjoys 
higher symmetries at the transition points. For instance, 
the XXZ model has a SU (2) symmetry at A = 1 while 
its symmetry group is weakened to SU q (2) for A ^ 1. 
However, when the transition takes place between an or- 
dered regime and a disordered one, the system remains 
gapless on the ordered side of the transition point and 
is gapful on the other side. Therefore, level-crossing of 
excited states is absent in this case. In our studies, the 
low- lying spectrum of the transverse-field Ising model has 
this character at A = 1. 

On the other hand, the concurrence measures actu- 
ally entanglement between two localized spins in a mixed 
state. Therefore, its behavior is also under the influ- 
ence of the excited states of the system, especially the 
low-lying excited state. To see this point more clearly, 
let us first consider a generalized model in one dimen- 
sion: H = J2i ( J x § i + JyS V i S i+i + J zSp z l+ i) + hsl 
which can be transformed into the ID XX Z model and 
the Ising model by appropriate choice of the parame- 
ters. We then introduce the identity: {0\[A, \H,A]]\0) = 
2J2 n ( E n - E o)\(0\A\n)\ 2 , where \n) and E n denote the 
eigenstate and the corresponding eigenvalue. The choice 
of operator A depends on the nature of local ordering. 
For the XX Z model near A = 1, antiferromagnetic or- 
der dominates so we take A = §" = J2j e y7r s" and obtain 

- £ W+i> = W + ivJ £o)|(0|.s»| 2 ,(4) 

a na 

where J = 2 + A. Obviously, the lhs of Eq. (gj is 
directly related to the concurrence of the XXZ model 
C = -2£ Q (s"s"+i) - 1/2. While for the Ising model, 
the ferromagnetic order dominates so A = §q = . s" 

an*) - w+i> - < j ^ + i> 

= ~ - ^j£(£« ~ £o)|<0|S»| 2 , (5) 

cm 

where J = — A. The lhs of Eq. JSJ) is also directly related 
to the concurrence of the Ising model C — 2((sfs|) — 
(sfsp — (sf§j)) — 1/2. The above two equations tell us 
that the concurrence does not simply depend on the re- 
scaled density of ground-state energy, but also the contri- 
butions from low-lying excited state with non-zero transi- 
tion amplitude of the order operator to the ground state. 

Keeping the above facts in mind, we now explain why 
the concurrence in the XXZ model reaches its maximum 
at A = 1. Notice that, on the left-hand side of this point, 
the lowest excited states are doubly degenerate and have 
spin numbers 5 = 1, S z = il. Correspondingly, the 
matrix elements of spin operators s+ and s~ between 
them and the singlet ground state |0), are nonzero. On 
the other hand, the second excited state has spin num- 
ber 5=1 and S z = and it contributes to the lon- 
gitudinal spin correlation function. On the right-hand 



side of the transition point, these excited states inter- 
change their position, as shown in Fig. [3 Since both 
sides are ordered phases, the main contributions to the 
corresponding order parameters are from the lowest ex- 
cited state. Then the rhs of Eq. J3J can be written as 
ffr + WjT, a (Ei- £ )|(0|J«|1)| 2 approximately. There- 
fore, the main contribution to the concurrence is from 
transverse order operator in the XY regime and from 
longitudinal order operator in Ising region. Only at the 
transition point A = 1, all three excited states have the 
same energy, and both the transverse and longitudinal 
spin correlation functions are power-law decay. Conse- 
quently, the contribution from both order operators to 
the concurrence makes it maximal. 

However, the above argument for the XXZ model is 
not valid for the transverse-field Ising model. For the lat- 
ter, if A > 1, the ground state is ferromagnetic and has 
parity P = 1 , while the lowest excited has parity P = — 1 . 
So the rhs of Eq. 10 is mostly contributed from the first 
excited state and the decreasing of the concurrence as A 
increases can be well understood. When A approaches the 
critical point, the gap formation in the thermodynamic 
limit will introduce a significant change to the concur- 
rence. This feature is reflected from the appearance of the 
minimum of the concurrence's first derivative at the crit- 
ical point. On the other hand, if A < 1, the phase is dis- 
ordered and gapful. In this situation, though the matrix 
element of order parameter between the ground state and 
the first excited state becomes smaller and smaller, the 
second excited state and other higher excited state now 
can not be neglected. Their participation not only com- 
pensates the lose from the first excited, but also makes 
the concurrence to be maximal at one point. However, 
when A — > 0, all excited state depart far away from the 
ground state which leads to the decrease of the concur- 
rence. Finally, for the Ising model, its singular behavior 
around the critical point just results from the transition 
from paramagnetic phase to ferromagnetic long-range or- 
der phase as discussed in Refs. [lLITcj. 

With this understanding, let us take another look at 
the Ji — Ji model. Besides the first-order transition point 
J2/J1 = 0.5, White and Affleck found a second-order 
transition at J%l J\ ~ 0.241 by numerical calculation [2l| . 
Around this point, the ground state of the model is non- 
degenerate but a level-crossing between two lowest ex- 
cited states occurs as shown in Fig. [21 Again, a similar 
equation of the concurrence involving matrix element of 
antiferromagnetic order parameter can be obtained. Un- 
like the XXZ model at A = 1, one of the excited states 
involved in level-crossing is a spin singlet and all the ma- 
trix elements of spin operators between it and the ground 
state, which is also a spin singlet, are zero. Consequently, 
this exicted state does not affect the concurrence at all 
so the concurrence does not show a maximum around 
J2/J1 ~ 0.241. Another excited state with spin 5=1 
involved in level-crossing is the only one who contributes 
to the concurrence significantly. Notice that this state is 
gapless for J2/J1 < 0.241 and is gapful otherwise. There- 
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fore, the transition is actually of order-to-disorder type, 
which we observed in the transverse-field Ising model. In 
this case we expect that the critical behavior of the con- 
currence will show up in its higher-order derivatives. To 
check this statement, we show minimum of higher order 
derivatives of the concurrence as functions of system size 
in Fig. [3J Although lack of sufficient data for careful 
finite size scaling analysis makes it unclear one of the 
minima will tend to J2 = 0.241, we believe that there ex- 
ists a certain order derivative of the concurrence whose 
minimum will tend to 0.241 for an infinite system. Such 
behavior is consistent with our picture. On the other 
hand, the concurrence is quite flat in the ordered phase 
and has its maximum at J2/J1 = (Fig. 1). This is 
due to the fact that the energy difference between the 
first excited state and the grounds state is almost a con- 
stant and the antiferromagnetic order parameter has its 
maximum at this point. 

Finally, for type II and III, we argue that the singu- 
larity of the concurrence at the critical point arises from 
the change of the long-range order. First, if the sys- 
tem does not have long-range order, then local prop- 
erties of the system, such as energy density, nearest- 
neighbor spin-spin correlation, etc., are not affected by 
spins far away, and their properties do not depend on 
the size of the system very much, nor does the concur- 
rence. On the other hand, if the system has long-range 
order, the correlation function in momentum space will 
be (5-function like, e.g., F(g)5(q—Tr) for antiferromagnetic 
long-range order, where F(g) is a coupling dependent 
function. Then the nearest-neighbor correlation function 



is, (0 1 erf a£ +s 1 0) cx F(g). Thus, the change of long-range 
order at the critical point, e.g., from paramagnetic phase 
to magnetic order phase for the Ising model Q], or from 
longitudinal magnetic order to that in the xy-plane for 
the 2D & 3D XXZ models will obviously result in the 
change of F(g) and leads to a singularity in the correla- 
tion function as well as the concurrence at the critical 
point. The size-dependent scaling behavior of the con- 
currence then comes naturally. 

In summary, based on the properties of the low-lying 
excitation spectrum near the quantum transition points, 
we examined and classified the critical behaviors of the 
concurrence, a measure of pairwise entanglement, in sev- 
eral typical spin models. We show that discontinuity of 
the concurrence in J\ — J2 model is simply caused by the 
ground-state level crossing at the transition point. On 
the other hand, QPT of the XXZ model at A = 1 is 
of ordered-to-ordered type and hence, is accompanied by 
a level-crossing between its lowest excited states. Conse- 
quently, its concurrence has a maximum at the transition 
point. Finally, for the transverse-field Ising model, whose 
transition at A = 1 is of ordered-to-disordered type, the 
opening of energy gap at the transition point introduces 
an extremum to the high-order derivative of the concur- 
rence. The singularity of the concurrence is the conse- 
quence of the transition of long-range-order. This classi- 
fication gives a unified picture on the connection between 
the QE and QPTs. 

This work is supported by RGC Projects CUHK 
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